Quantum Confinement Induced Molecular Mott Insulating State in La4Ni 3 Os 
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The recently synthesized layered nickelate La4Ni 3 0s, with its cuprate-like Ni02 layers, seemingly 
requires a Nil (d 8 )+2Ni2 (d 9 ) charge order, together with strong correlation effects, to account for 
its insulating behavior. Using density functional methods including strong intra-atomic repulsion 
(Hubbard U), we obtain an insulating state via a new mechanism: without charge order, Mott 
insulating behavior arises based on quantum coupled, spin-aligned Ni2-Nil-Ni2 d z i states across the 
trilayer (rather than based on atomic states), with antiferromagnetic ordering within layers. The 
weak and frustrated magnetic coupling between cells may account for the small spin entropy that 
is removed at the Neel transition at 105 K and the lack of any diffraction peak at the Neel point. 



Nearly 25 years since the discovery of high tem- 
perature superconductivity (HTS), the question per- 
sists: why does HTS occur only in copper oxides and 
not other layered oxides. Nickelates with cuprate- 
like structures [1[ have seemed to be the most likely 
candidates: the Cu 2+ d 9 configuration with active 
d x 2_ y 2 orbital maps onto the Ni 1+ configuration, 
but this is known to be a problematic charge state. 
Nevertheless, difficult charge states may sometimes 
be stabilized, such as the nominal Ni 1+ ("infinite 
layer") compound LaNi02 prepared by topotactic 
exchange of oxygen Q. Cuprate-like electronic struc- 
tures in other oxides have been pursued since soon 
after HTS were discovered. The charge conjugate 
analog S^VCU (d 1 vs. d 9 ) presented hopel3|, |4[ 
and displays magnetic and orbital ordering [5] but 
has never shown superconductivity. Layered nicke- 
lates have been seen as the best mimic of HTS if the 
d 9 Ni + charge state can be stabilized, and theoret- 
ical interpretation of the few known prospects has 
begun@,0|. The LaNi0 3 /LaA10 3 superlatticefl, ©| 
is an example that has been suggested to produce 
an HTS-like Fermi surface, but is based on the more 
common Ni 3+ ion and has not yet produced super- 
conductivity. 

The discovery by Greenblatt's group of a 
Ruddleston-Popper sequence [10] of phases 
La n+ iNi n 02n+2 with n cuprate-like MO2 lay- 
ers has reinvigorated interest in nickelates, with a 
focus on L^NiflOg ("La438") where stoichiometry 
is attained [111. Il2j . In an ionic picture, the mean 
Ni valence in La n+ iNi n 02 n +2 is +(n + l)/n, i.e. 
metallic. However, La438 is highly insulating [131] . 
with a room temperature resistivity ~ 10 3 O-cm 
increasing by six orders of magnitude down to 25 
K. Magnetic ordering is observed[13j at Tat=105 
K, with NMR data indicating an antiferromagnetic 
(AFM) transition^. The entropy loss at Tat 
is only ~Rln2 per Ni, suggesting unconventional 
magnetic behavior in this system. 



La438 crystallizes in a tetragonal unit cell, with 
space group I4/mmm (#139) and lattice parameters 
a=3.9633 A c=26.0373 A. The structure, pictured 
in Fig. [TJ consists of three MO2 "infinite layer" 
planes separated by layers of La 3+ cations but no 
oxygen. On either side of this trilayer slab lies a flu- 
orite La/02/La layer, with the same structure as the 
corresponding layer in LaFeAsO. While in a sense 
isoelectronic with doped cuprates as it seems to in- 
volve some mix of Ni 1+ and Ni 2+ cations, the elec- 
tronic and magnetic structure of the Ni ions is still 
unclear. As mentioned, the low oxidation state Ni 1+ 
is enigmatic; to obtain such a configuration for Ni 
is not easy, and these Ruddleston-Popper phases are 
unstable at moderate (~375 °C) temperature [lOj to 
decomposition into La203 and Ni metal. In the Ni 1+ 
"infinite-layer" compound LaNi02 it was found to 
be impossible [7] to produce a Mott insulating state 
as might be anticipated; in fact the compound is 
metallic at least roughly as calculated. La438 con- 
sists of three layers of LaNi02 subject to quantum 
confinement, as we demonstrate below. 

In La438 a simple ionic picture will give Ni con- 
figurations d 8 + 2d 9 , which would give the natural 
identification of the one 'inner' Nil = Ni 2+ (d 8 ) and 
two outer Ni2 = Ni 1+ (d 9 ) cations, as labeled in Fig. 
[TJ Nil-Ni2 layers are separated by only 3.25 A, much 
less than the in-plane separation a=3.93 A, a char- 
acteristic that will become important below. The 
Ni trilayers are separated by the La/02/La unit, 
and an (a/2, a/2,0) translation between successive 
MO2 trilayers reinforces very small interplanar cou- 
pling between cells, which will be shown to produce 
a natural quantum confinement within the trilayers. 

Here we study the electronic structure of the com- 
pound by ab initio techniques to understand the 
electronic and magnetic structure of La438 and to 
identify the differences and similarities with respect 
to the cuprates. Our electronic structure calcula- 
tions were performed within density functional the- 



ory [15[ using the all-electron, full potential code 
wien2k [l6| based on the augmented plane wave 
plus local orbital (APW+lo) basis set [17], and 
the experimental atomic positions [TjJ. The gener- 
alized gradient approximation [18] was used. To deal 
with str ong correlation effects we apply the LDA+U 
scheme [jjl EqI including an on-site repulsion U and 
Hund's coupling J for the Ni 3d states. Results 
were not very dependent on the specific values of 
U and J in the reasonable range, and we report re- 
sults with U = 4.75 eV, J = 0.68 eV, values very 
similar to those determined from constrained den- 
sity functional calculations fl3j| . 

The bands within LDA (U=0) are metallic as ex- 
pected, and in fact the bands crossing the Fermi 
energy Ep are dispersive (bandwidth of ~2-3 eV). 
These bands do not suggest substantial effects of 
strong correlation, and Poltavets et al. sugg ested[13] 
Fermi surface driven spin density wave (SDW) order, 
which might destroy part of the Fermi surface and 
account for the magnetic transition. The earlier cal- 
culations, even those using LDA+U, produced only 
metallic states [13]. The key to obtaining a lower 
energy, and insulating, state is to allow extra free- 
dom by doubling the lateral cell of La438, and we 
found that in-plane AFM order is readily obtained. 
Various configurations and types of spin order were 
tried, but no inter-plane AFM order could be ob- 
tained; moments along the z axis always want to 
align within the trilayer. 

We discuss only this most energetically stable 
state, comprised of AFM-ordered Ni layers with 
spins aligned along the z axis. The spin density, 
pictured in Fig. [TJ makes evident both the type of 
spin ordering and the character (shape) of the spin 
density, which reveals both d x i_ y i and d z i contribu- 
tions. There is indeed some charge difference with 
muffin-tin charges of 8.52 (Nil) and 8.64 (Ni2), and 
different moments of 1.4 and 1.3 fiB respectively. 
It is understood that different formal charge states 
may involve surprisingly reduced difference in ac- 
tual charge; it is magnetic states that are more in- 
structive. While the moments we calculate might be 
consistent with an S=l configuration (d 8 ) reduced 
by hybridization, neither moment (1.3-1.4 (ib) can 
be ascribed to an S=^ d 9 ion (which would have 
a maximum moment of 1 /i#). Without d 9 ions, 
the ionic, charge-ordered picture becomes untenable, 
and a new mechanism for producing the gap must 
operate. 

Our LDA+U band structure, shown in Fig. [2j pro- 
duces a Mott insulating state, but of an unusual type 
that requires broadening one's view of how these cor- 
relation effects operate. In the AFM state, the spin- 




FIG. 1: Left side: structure of I^NisOs, showing the 
three NiCb "infinite layer" structures sandwiched on ei- 
ther side by the fluorite La/02/La blocking layer. The 
Nil and Ni2 designations of the NiCb trilayer are pro- 
vided. Right side: spin density isocontour plot of the 
ground state of La4Ni30s from LDA+U calculations for 
the a/2 x a/2 x 1 cell. The different colors denote the two 
different spin directions. Magnetic order is AFM within 
layers, with spins aligned along the perpendicular (ver- 
tical) axis. 



up and -down band structures are of course identical. 
Since the valences of La 3+ and O 2- are clear, the Ni 
ions are on average Ni 1,33+ , with charge d 8 67 mak- 
ing the maximum moment around 1.33 which is 
in fact what they are in our ground state. The Ni 
ions are thus fully polarized, meaning the majority 
orbitals are completely occupied and removed from 
consideration. Also, the minority d x i_ y i orbitals are 
empty, centered at +3.5-4 eV on the energy scale of 
Fig. [2j Thus the d x i_ y i orbital contributes one unit 
to the hole count, 1 \±b to each Ni moment, and the 
in-plane superexchange coupling, but otherwise is of 
no interest. The relevant orbitals are solely the par- 
tially occupied minority d z i orbitals on the six Ni 
ions in the AFM cell, which must accommodate one 
hole per triple (two in the AFM cell). The central 
question then is: how can a Mott insulating state, 
as we obtain, arise at "2/3 filling" if charge dispro- 
portionation into 2 Ni 1+ + Ni 2+ does not occur? 

The in-plane AFM order, as is common, drasti- 
cally reduces bandwidths (spin-conserving hopping 
takes place only between second Ni neighbors in 
the plane), so taking account of Ni d interplanar 
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FIG. 2: The Mott insulating band structure of La4Ni3 08 for the calculated ground state of Fig. U Left panel: the 
basic plot, allowing nearly degenerate bands to be distinguished. Other panels: I^NisOs fatband representation 
emphasizing (middle panel) Ni2, and (right panel) Nil, d character. Note that the band at the bottom of the gap 
has vanishing Nil character. 
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FIG. 3: Energy level diagram. The intraatomic levels are shown on the left, with occupied majority states on the 
far left, partially occupied minority levels just to the right. The right panel indicates how the d z i orbit als on the 
Ni2-Nil-Ni2 triple couple to give a gap at 2/3 band filling, following Eq. [1] 



coupling t± (between Nil and Ni2) becomes impor- 
tant. For each Ni2-Nil-Ni2 triple along the z axis 
with neighbors coupled by t± (with zero site energy, 
neglecting the small on-site energy difference), the 
eigenvalues and eigenvectors are 

Ej = 0, ±y/2t±, 
\E 3 > = -1(1,0,-1), 1(1,±V2,1). (1) 

The non-bonding (odd symmetry) state at Ej = q=0 
which does not involve the Nil site is flanked be- 
low and above by even symmetry bonding (E+) and 
antibonding (E-) combinations. This coupling is 
pictured at the right of Fig. [3] 

In the middle and right panels of Fig. [2] the Ni 
d characters are emphasized for Ni2 and Nil, re- 
spectively. The bands display both the spectrum 



(most clearly evident near the M point) and the 
(Nil vs. Ni2) character of this "triatomic molecule" . 
For example, the band at the bottom of the gap is 
(non-bonding) purely Ni2, which is Eq, while the 
other two bands have both Nil and Ni2 character 
(E+ , E_ ) . Two of these levels will be filled (below 
the Fermi level) and one empty (above the Fermi 
level). The gap size, which is 0.6 eV here, will de- 
pend both on the coupling t± and the value of U. 

Since previous LDA+U studies produced only a 
metallic result, one can understand the electronic 
structure only using a quantum-coupled unit, a 
molecular trimer basis, for the trilayer, from which 
on-site repulsion produces a (Mott) insulating state 
in agreement with observation [T^j. The intralayer 
AFM magnetic coupling is strong and is mediated 
by d x 2_ y 2 superexchange, as in cuprates. Between 



Ni02 layers, the coupling is FM direct exchange aris- 
ing from the d z 2 orbitals (Wannier functions). This 
FM coupling is robust in the sense that we could 
not obtain antialigned layer magnetic order. Exper- 
imental data suggest an in-plane magnetic coupling 
with J || « 120 K[13], which is substantial but an 
order of magnitude smaller than in cuprates. 

The trilayer-to-trilayer coupling will be superex- 
change through the fruorite La/02 /La spacer layer 
that can be seen in Fig. [TJ The k z bandwidth and 
related dispersion can be estimated from the split- 
ting of the band below the gap near -0.5 eV at T (this 
splitting vanishes at tt/c). The hopping t z ~ 15 meV 
from a Ni2 ion to four Ni2 ions across the blocking 
layer leads to J z = U 2 JU ~ 0.2 meV - 2 K. Such 
small interlayer coupling enters T^r logarithmically 
and can still give rise to magnetic order. However, 
the interslab coupling is frustrated as well as being 
small. These facts suggest that some of the entropy 
will not be lost at Ttv=105 K (full magnetic order 
along the z axis is not established), consistent with 
the small observed entropy change and also with the 
absence of a new diffraction peak[l3| below Tat. 

The state we have obtained and analyzed has in- 
sulating, AFM order in-plane in common with the 
HTS cuprates, as well as the obvious similarities in 
structure. The differences are however substantial. 
First, excellent insulating character and relatively 
strong superexchange-mediated AFM order is ob- 
tained without integer formal oxidation states on the 
Ni ions. Such behavior has not been reported previ- 
ously for transition metal oxides. Second, the d x i_ y i 
hole on each Ni ion is robust, and doped carriers go 
only into d z i states; hence the AFM order (driven by 
d x 2_ y 2 coupling) may survive to much larger doping 
level than in cuprates. Orbital order is ruled out. 
Third, holes will enter the non-bonding band that 
is confined to the outer Ni2 layers, so beyond the 
insulator-to-metal transition these two conducting 
layers will be separated by a non-participating Nil 
middle layer; the conductivity will be both orbital 
and layer specific. Doped electrons will go into all 
three Ni layers, but with half of the charge on the 
central Nil layer. 

The insulator-to-metal transition in either case 
will likely occur within the AFM ordered phase 
rather than occurring simultaneously. Superconduc- 
tivity coexisting with AFM order is known to occur, 
but the competition tends to be very unfavorable 
for superconductivity. Moreover, the magnetic in- 
teractions between the doped carriers (in d z i states) 
will be very weak, and very unlike cuprates where 
carriers are doped into the d x i_ y i states. Though 
substantially different from cuprates, doped La438 



should be an extremely interesting system, and dop- 
ing should be pursued. 

We acknowledge informative communication with 
N. J. Curro. This project was supported by DOE 
grant DE-FG02-04ER46111. 



* Electronic address: victor. pardo@us cTes] 

* Electronic address: wepickett@ucdavis.edu 

[1] W. E. Pickett, Rev. Mod. Phys. 61, 433 (1989). 
[2] M. Kawai, S. Inoue, M. Mizumake, N. Kawamura, 

N. Ichikawa, and Y. Shimakawa, Appl. Phys. Lett. 

94, 082102 (2009). 
[3] M. Cyrot, B. Lambert- Andron, J. L. Soubeyrous, 

M. J. Rey, P. H. Dehauht, F. Cyrot- Lackmann, 

G. Fourcaudot, J. Beille, and J. L. Tholence, J. Solid 
State Chem. 85, 321 (1990). 

[4] D. J. Singh, D. A. Papaconstantopoulos, 

H. Krakauer, B. M. Klein, and W. E. Pickett, 
Physica C 175, 329 (1991). 

[5] H. D. Zhou, B. S. Conner, L. Balicas, and C. R. 

Wiebe, Phys. Rev. Lett. 99, 136403 (2007). 
[6] V. I. Anisimov, D. Bukhvalov, and T. M. Rice, Phys. 

Rev. B 59, 7901 (1999). 
[7] K.-W. Lee and W. E. Pickett, Phys. Rev. B 70, 

165109 (2004). 
[8] J. Chaloupka and G. Khaliullin, Phys. Rev. Lett. 

100, 016404 (2008). 
[9] P. Hansmann, X. Yang, A. Toschi, G. Khaliullin, 

O. K. Andersen, and K. Held, Phys. Rev. Lett. 103, 

016401 (2009). 
[10] V. V. Poltavets, K. A. Loshkin, T. Egami, and 

M. Greenblatt, Mat. Res. Bull. 41, 955 (2006). 
[11] V. V. Poltavets, K. A. Loshkin, M. Croft, T. K. 

Mandal, T. Egami, and M. Greenblatt, Inorg. 

Chem. 46, 10887 (2007). 
[12] V. V. Poltavets, M. Greenblatt, G. H. Fecher, and 

C. Felser, Phys. Rev. Lett. 102, 046405 (2009). 
[13] V. V. Poltavets, K. A. Loshkin, A. H. Nevidomskyy, 

M. Croft, T. A. Tyson, J. Hadermann, G. V. Ten- 

deloo, T. Egami, G. Kotliar, N. ApRoberts- Warren, 

et al., Phys. Rev. Lett. 104, 206403 (2010). 
[14] N. ApRoberts- Warren, A. Dioguardi, V. V. 

Poltavets, M. Greenblatt, and N. J. Curro, 

arxiv/1005.0435 (2010). 
[15] R. O. Jones and O. Gunnarsson, Rev. Mod. Phys. 

61, 689 (1989). 
[16] K. Schwarz and P. Blaha, Comp. Mat. Sci. 28, 259 

(2003). 

[17] E. Sjostedt, L. Nordstrom, and D. J. Singh, Solid 

State Commun. 114, 15 (2000). 
[18] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. 

Rev. Lett. 77, 3865 (1996). 
[19] V. I. Anisimov, J. Zaanen, and O. K. Andersen, 

Phys. Rev. B 44, 943 (1991). 
[20] E. R. Ylvisaker, W. E. Pickett, and K. Koepernik, 

Phys. Rev. B 79, 035103 (2009). 



4 



